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Abstract 

The notion of differentiation index for DAE systems of arbitrary order with generic 
second members is discussed by means of the study of the behavior of the ranks of 
certain Jacobian associated sub-matrices. As a by-product, we obtain upper bounds 
for the regularity of the Hilbert-Kolchin function and the order of the ideal associated 
to the DAE systems under consideration, not depending on characteristic sets. Some 
quantitative and algorithmic results concerning differential transcendence bases and 
induced equivalent explicit ODE systems are also established. 
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1 Introduction 



This paper is devoted to the study, mainly from a quantitative point of view, of differential 
algebraic equation (DAE) systems of the form: 

h(X,U) = Xt 

(y\ ._ J fn(X, U) = X n 

1 ' ' I gi(X,U,U,...,U^) = Y X ' 

, g r (X,U,U,...,U^)) = Yr 

where fi, ■ ■ ■ , f n are polynomials in the n+m variables X := X\, . . . , X n , U := U\, . . . , U m , 
and, for every 1 < j < r, gj is a polynomial in the n + (ej + l)m variables X, U and 
the derivatives U"> := U{ , . . . , U m , 1 < i < Bj, with coefficients in a differential field k 
(for instance := Q, M, C, Q(t), etc.). The constants ej S No denote the order of the 
respective equation gj in the variables U. The variables Y := Y\, . . . , Y r are a new set of 
differential indeterminates which can be viewed as parameters (while the variables X and 
U are the unknowns of the system). So, it is quite natural to extend the ground field k to 
the differential field L := k(Y) (i.e. the smallest field containing k and all the successive 
derivatives of Y) and consider our input system also as a system over L. Even if we do not 
assume, as customarily, a differential O-dimensional situation (r may be strictly smaller 
than m), we will suppose that the last r equations are "independent" in a suitable natural 
way defined in Section 12.11 

DAE systems like (£) can be regarded from several points of view: for instance, this 
kind of systems arises in Control Theory (see for instance [8] and [3l Section 4]); they may 
also be interpreted as the equations defining the graph of a differential morphism (see 
|21|). The system (£) may be viewed as a family of usual polynomial DAE systems where 
the second member parametrizes the family and takes arbitrary values outside a suitable 
proper algebraic Zariski closed set (see (30j Section 5.2]). In this last sense we say that 
the system (£) is generic. 

We will focus on several basic topics concerning the system (£), including the number 
of differentiations that suffice to obtain explicit equations, a description of all the relations 
of a prescribed order that all the solutions must verify and the number of initial conditions 
that can be arbitrarily fixed. All these aspects have been studied extensively during the 
last two centuries beginning with two posthumous articles by Jacobi ( [12] and [13]). The 
present paper (as well as most of the previous ones on the subject) may be considered as 
a modern approach to the work done by Jacobi in these remarkable and not sufficiently 
known papers. 

The main notion we will consider is the differentiation index, which is a well known 
and important invariant associated to a DAE system (customarily for first-order and 0- 
dimensional systems). 

There are many different, not always equivalent, definitions of differentiation indices 
(see for instance [§], [H], [3], [5], [27], [2S], [TB], [2U], [33], [23] and [21]). Here we are 
mainly interested in the so-called global differentiation index (see [21 Section 2.2]). Roughly 
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speaking, the differentiation index represents the minimum number of times that all or 
part of a given DAE system must be differentiated in order to obtain an equivalent explicit 
ordinary differential equation ( ODE) system; in other words, the number of differentiations 
required to determine the derivatives of a certain order of the unknowns as continuous 
functions of derivatives of lower order of the unknowns (see [21 Definition 2.2.2]). In some 
sense, the index can be regarded as a measure of the complexity of the DAE system. 
From the theoretical point of view it represents the distance between the given system 
and another one for which an existence and uniqueness theorem holds (see |26j). On the 
other hand, from the point of view of its numerical resolution, it is closely related to the 
condition number of the iteration matrix in the implicit Runge-Kutta method (see [21 
Theorem 5.4.1]). 

In this paper we give a precise algebraic definition of a differentiation index for DAE 
systems as (£), not necessarily of first order nor ^-dimensional (see Definition [TBI below) . 
by means of certain stationary properties of the ranks of suitable Jacobian sub-matrices 
which are proved in Section [3l Another equivalent definition of this index, in terms of 
a quite natural filtration given by the successive differentiation of the input equations, 
is implicitly contained in Theorem [T7] below. In particular, this last formulation shows 
how the differentiation index uncovers constraints that every solution must satisfy. This 
approach is closely related to the classical algorithmic definitions of the differentiation 
index by means of iterated prolongations (differentiations) and projections (eliminations) 
(see [26], 0, [IS] and [21]). 

Our theoretical approach has the advantage of leading to a polynomial time algorithm 
which computes the index of the input system (£) by simple comparison of ranks of 
Jacobian matrices (see Section I5.3.2P . For previous work related to the computation of 
differentiation indices of DAE systems, we refer the reader to [5], [18], [26], [21], [25], [23] , 
[21], P5 and [J9]. 

As we said before, the notion of differentiation index is closely related to the possibility 
of writing certain derivatives of the unknowns in terms of derivatives of lower order of these 
unknowns. Unfortunately, in general, this cannot be done using only the algebraic relations 
induced by the original equations. For example, if all the equations have order one, this 
nice situation, which is usually ensured by the Theorem of Implicit Functions, corresponds 
exactly to those systems whose associated index is 0. However, in the general case, by 
successive differentiations (as many as the index) we can always obtain such a situation. 
Evidently, the new explicit system comes out of the frame of the polynomial (or even 
rational) systems. Under certain additional conditions on the system (£), we are able to 
compute an implicit ODE system which is equivalent to (£) (see Section[5]); in particular, in 
the case of a first-order system, we give an implicit but simple polynomial way to describe 
it, distinguishing the variables by their interrelations (namely, "free variables", "implicit 
variables" , etc.). Moreover, we can estimate degree and order upper bounds for the implicit 
equations involved in the equivalent ODE systems (see Subsection I5.3.1f) and give an 
algorithm to compute them with polynomial complexity in terms of an intrinsic parameter 
related to the geometric degree of suitable associated algebraic varieties (Proposition l33l) . 

Our method also allows us to give a new upper bound for the regularity of the Hilbert- 
Kolchin function (or differential Hilbert function) associated to the DAE system (£). The 
Hilbert-Kolchin function is introduced in [141 Chapter II] in order to estimate, for each 
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non- negative integer i, the degree of freeness of the first z-derivatives of the unknowns 
modulo the relations induced by the input equation system (see also Section 1131 below for 
a precise definition in our case). 

As it happens for the classical Hilbert function associated to homogeneous polynomial 
ideals (see for instance [TJ Chapter 11]), the Hilbert-Kolchin function becomes a well 
defined polynomial for sufficiently big arguments i (in the ordinary differential setting, 
this polynomial is extremely simple since its degree is at most 1). The regularity of the 
Hilbert-Kolchin function is defined to be the first non-negative integer from where the 
function and the polynomial coincide. It is well known that this regularity can be exactly 
described in terms of the orders of the elements in a characteristic set associated to any 
orderly ranking (see the proof of [141 Chapter II, Section 12, Theorem 6] or [U Theorem 
3.3]). In this article (see Theorem 1191 below) we show that max{l, e%, . . . , e r } — 1 is an 
upper bound for the regularity of the Hilbert function for the system (£) over the field 
L (notice that no characteristic sets are involved in the bound). In particular, if (£) is a 
first-order system, the regularity of the Hilbert-Kolchin function is 0; in other words, the 
Hilbert-Kolchin function and the associated polynomial coincide for all i. 

As further consequences of our techniques we deduce Bezout-type bounds for the differ- 
entiation index and the order of the input system (X) in terms of the orders of its defining 
equations (see Remark 1161 and Subsection 14.2.2] ). A similar bound for the differentiation 
index may also be obtained by rewriting methods as in [30] Section 5.2] (see also Section 
3] for an analogous bound for 0-dimensional first-order systems). Concerning the order, 
we recover Ritt's differential analogue to Bezout's Theorem ([28] Ch. VII, p. 135]). A more 
precise bound for the order of a 0-dimensional system was conjectured by Jacobi (see |12j ) 
and proved in [15] under additional hypotheses which are met in our situation (see also 
[29] . [22]). We point out that a convenient refinement of our approach may be applied to 
obtain Jacobi-type bounds for both the order and the differentiation index, which will be 
the subject of a forthcoming paper. 

The paper is organized as follows: basic definitions and notations are introduced in 
Section [2j Section [3] is devoted to the study of the behavior of an integer sequence which 
is strongly related to the ranks of suitable Jacobian sub-matrices of the input equations 
and their successive derivatives. In Section H] we give a precise algebraic definition of the 
differentiation index based on the behavior of this sequence, and an equivalent description 
of it in terms of the variety of constraints (Subsection 14. ip . Our results on the Hilbert- 
Kolchin function of the differential ideals associated with the considered DAE systems are 
presented in Subsection 14. 2t we estimate the regularity of the Hilbert-Kolchin function in 
Subsection 14.2.11 and the order of the induced differential ideal in terms of the orders of the 
equations in the given system in Subsection l4.2.21 further, a first result related to equivalent 
explicit ODE systems is given for the differential 0-dimensional case in Subsection 14.2.31 
This result is refined in Section [5] by introducing a more accurate definition of the index, 
which is done in Subsection 15.1] and generalized for first-order systems with positive 
differential dimension in Subsection 15.21 Finally, we present quantitative and algorithmic 
considerations concerning these results in Subsection 15.31 

Acknowledgements The authors thank Evelyne Hubert (INRIA, Sophia Antipolis), 
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Gustavo Massaccesi (Universidad de Buenos Aires), Francois Ollivier (Ecole Polytech- 
nique, Palaiseau) and Alexandre Sedoglavic (Universite de Lille) for their helpful remarks. 

2 Preliminaries 

2.1 Basic definitions and notations 

Let A; be a characteristic zero field equipped with a derivation S (for instance k = Q, R or 
C with 5 := 0, or k = Q(t) with the usual derivation S(t) = 1, etc.). 

For an arbitrary set of (differential) indeterminates Z := Z\, . . . , Z a over k we denote 
the Z-th successive derivative of a variable Zj as Z^ (as customarily, the first derivatives 

are also denoted by Zj); we write Z® := {zf , . . . , Z^} and Z® := {ZW, < i < I}. The 
(non-noetherian) polynomial ring k[Z^ l \ I £ No], called the ring of differential polynomials, 
is denoted by k{Zi, . . . , Z a } and its fraction field by k(Z\, . . . , Z a ). Given a finite set of 
(differential) polynomials H±, . . . , Ha € k{Z±, . . . , Z a }, we write [Hi, . . . , Hp] to denote 
the smallest ideal of k{Z\, . . . , Z a } stable under differentiation, i.e. the smallest ideal 
containing Hi, . . . , Hp and all their derivatives of arbitrary order. The ideal [Hi, . . . , Hp] 
is called the differential ideal generated by Hi, ... , Hp . 

We deal with a particular class of differential algebraic equation (DAE) systems: 

fi(X,U) = Xi 

fn{X,U) = X n . . 

gi(X, U,U,..., U^) = Yi { ' 

, g r {X,U,U,...,U^)) = Yr 

where fi, ■ ■ ■ ,f n are polynomials in the n+m variables X := X\, . . . , X n , U := U\, . . . , U m , 
and, for every 1 < j < r, gj is a polynomial in the n + (ej + l)m variables X, U and 
the derivatives U (l> := U{ , . . . ,Um , 1 < i < ej, with coefficients in the field k. The 
constants ej € No denote the order of the respective equation gj in the variables U, i.e., 
the order of the highest derivative of a variable in U appearing with non-zero coefficient 
in this polynomial. The variables Y := Yi, . . . , Y r form a new set of indeterminates which 
we regard as parameters (while the variables X and U are the unknowns of the system). 
We allow n to be equal to zero (i.e. no variables X appear in the system). 

In addition, we will assume that the polynomials gi are differentially algebraically 
independent over k as elements of the fraction field Frac(/c{y, X, U}/[f\ —X\, . . . , f n —X n ]), 
i.e. there is no non-trivial differential relation involving the classes of the differential 
polynomials gi,...,g r over k. This assumption guarantees that the differential ideal of 
k{Y, X,U} associated to system ([TJ does not contain a non-zero differential polynomial 
involving only variables Y. 

For any differential polynomial g lying in a differential polynomial ring k{Zi, . . . , Z a } 
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the following recursive relations hold for the successive derivatives of g: 
g® := 5{g^)) + Y, ij ^^zf +1 \ for I > 1, 

,3 gz (l) J 

where 5(g^ l ~ 1 ^) denotes the polynomial obtained from g( l ~ 1 ^ by applying the derivative 5 
to all its coefficients (for instance, if A; is a field of constants, this term is always zero). 
Concerning the system ([T]) we also introduce the following definitions and notations: 

• L denotes the fraction field k(Y). 

• e := max{l, ei, e2, • • • , e r }. 

• For every I G No we set: 

if) := /(0 G k[xWu®] i = l,...,n, 
G f := gf - Y® G k[Y M , X® , LT^] j 1 r. 

• For every Z G N, A\ denotes the polynomial ring A\ := L[xW, U"*] and A^ C A/_x+ e 
the ideal generated by ^, G^ 1 ! (observe that the ideal A/ is contained in Ai_i +e 
because the orders of i^ -1 ) and G^" 1 ) are at most Z — 1 + e). We set Ao := (0). 

• A := [F, G] C U} is the differential ideal generated by the polynomials F := 
F u ...,F n andG:=Gi,...,G r . 

• IK is the differential field k(X)(U) with the derivation induced by Xj = fi(X,U), 
i = 1, . . . , n. 

The hypothesis on the differential algebraic independence of the polynomials gj (1 < 
j < r) easily implies that t/ie ideaZ A has differential dimension m — r over L (see for 
instance Proposition 12]). Furthermore, the map 

Xj ^ fj, J = 1, • ■ ■ , n, Yj ^ gj, j = 1, . . . , r, and ?7j h-> C^-, j = 1, . . . , m, 

induces an isomorphism between the fraction field Frac (L{X, U}/A) and the differential 
field K = k(X)(U) (see Remark 7]). 

In the sequel, in order to simplify notations, for any g G L{X, U}, we will also write g 
for its class in L{X, U}/A or in its fraction field IK (the ring where we consider the object 
will be clear by the context). In the same way, g will denote either the derivative of g in 
L{X, U} or its derivative as an element of IK, and so on. 

2.2 Associated Jacobian sub-matrices 

Here, we introduce the Jacobian matrices and sub-matrices we will deal with throughout 
the paper. From the considered differential system ([TJ), we define a family of sub-matrices 
constructed from the (infinite) Jacobian matrix associated to the (infinitely many) poly- 
nomials F® and G® with respect to the (infinitely many) variables and U^>. 
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Definition 1 For each k £ N and i £ N> e „i (i.e. j 6 Z and i > e — 1), we define the 
matrix Zk,i G K k( - n+r ^ xk( - n+m ^ as follows: " 

9t/(i+l) U U U U 

g G (i- e +l) ggji-e+l) 



Qp(i-e + k) Qp(i-e + k) Qp(i-e+k) Qp(i~e+k) Qp{i-e + h) Qp(i-e+k) 

QQ(i-e+k) QQ(i-e + k) Q(j(i-e+k) QQ(i-e + k) Q(j(i-e+k) QQ^-e+k) 

\ ax(*+i) ~dUU+ T T ~dxTi+W a(/(»+2) ' ' ' dX( i + k ) / 
In other words, 3fc,i is the Jacobian matrix of the polynomials 

p{i— e+l) gr(i— e+l) p(i—e+k) Q(i—e+k) g e+fe+1] jy[«+fc]j 

with respect to the variables C/^ +1 ^, . . . where the entries are re- 

garded as elements in K. 

Observe that the block triangular form of Zki follows from the fact that the differential 
polynomials i^(*~ e +0 and G^~ e+l ^ have order bounded by i + I. Hence, their derivatives 
with respect to the variables and are identically zero for j > I + 1. 

The matrices Zk,i are strongly related with some algebraic facts concerning the (alge- 
braic) ideals A/ introduced in the previous subsection: 

Proposition 2 Let k £ N and i £ N> e _i. Then: 

(i) The transcendence degree of the field extension 

Frac(Ai/(Ai_ e+1+fc n At)) ^ Fvac{A i+k / Aj_ e+1+fe ) 

equals the dimension of the kernel ofZk,i- 

(ii) The following identity holds: 

trdeg L (Prac(Aj/Aj_ e+1+fc n Ai)) = (m - r){i + 1) + e(n + r) - dim K (ker(^ J), 
where Z\ j denotes the usual transpose of the matrix Zk,i- 



Proof. In order to prove (i) we follow closely the proof of Proposition 16] (see also 
Theorem 10]): since the polynomials and Gj have order p+1 and p+ej respectively, 
we conclude that for i > e — 1, the polynomials F^ l ~ e ^ and G^~ e ^ belong to the ring A{. In 
any case (even if i = e — 1), we may consider Frac(j4j + fc/Aj_ e+1+ / c ) as the fraction field of 
the domain 

R := K[X^ i+l \ U^ i+l \ . . . , X^ i+h \ J/( i+fc )]/(i? 1 ( ,; ~ e + 1 ) ) G^~ e+l \ . . . , F^~ e+k \ G^~ e+k h, 
where if denotes the field Frac(^4j/(Aj_ e+ i + fc n Ai)). 
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Then, the transcendence degree we want to compute is the dimension over the field 
Frac(-R) (or equivalently, over Fiac(Ai + k/ ^i-e+i+k)) of the kernel of the Jacobian matrix 
associated to the K— algebra R (see for instance \\.7\ Chapter VI, §1, Theorem 1.15]). This 
Jacobian matrix is exactly the matrix Zk,i- To finish the proof of (i), it suffices to show 
that the dimension of both kernels (namely, over Frac(i?) and over K) of Zk,i is the same. 

Let us observe that the entries of Zk,i can De regarded as polynomials in the ring 
k[X, U^} C K, which is isomorphic to k[Y^~ e+k \ X^+ k \ U^]/(F^ e+k \ G^~ e+k ^). 
After tensoring by L, this last ring becomes Ai +k / A-i-e+i+k- Now, since LnA = 0, the rank 
of Zk,i is preserved after that tensoring and so, its rank over the fraction field of k[X, U^} 
(and therefore, over K) is equal to its rank over the fraction field of Ai +k / £±i- e +i+k 
(namely, Frac(i?)). This finishes the proof of (i). 

For part (ii), it is easy to see that the polynomials F^ P \G^ form a regular sequence 
(see for instance Corollary 9]) and hence, the field extension L <— > Frac(^4j + fc/Aj_ e _|_i + fc) 
has transcendence degree equal to the number of variables in Ai +k minus the number of 
generators of the ideal Aj„ e+1+/ t, that is (i + k + l)(n + m) — (i — e + k + l)(n + r) = 
(m — r)(i + k + 1) + e(n + r). The result follows by considering the tower of fields 

L ^ Frac(Aj/Ai_ e+ i +fc n At) ^ Fvac(A i+k / Ai^ e+1+k ) 

and part (i), noticing that diniK kei(Zk,i) = k(m — r) + dimuker^ f ). ■ 

Remark 3 The equality in part (ii) of Proposition [H also holds for k = defining 
diniK(ker(3o A) := for i > e — 1. This follows from the fact that Aj_ e+ i C A{ is 
generated by a regular sequence, as in the proof of that Proposition. 

Remark 4 Fix an index i > e — 1. Note that the prime ideals (Aj_ e +i+fc n Ai)k^N form 
an increasing chain of ideals. Since Ai is a noetherian ring, this chain must be stationary 
for k big enough, and so, from condition (ii) of Proposition^ the sequence of integers 
(dims ker(3t j))fceN ^ s non- decreasing and becomes stationary for k big enough. 

In the next section we proceed to study more closely this kind of stationary properties 
related to the rank of the matrices Zk,i- 

3 The rank of Jacobian sub-matrices 

The aim of this section is to study the behavior of the ranks of the matrices Zk i when 
k and/or i run over N and N> e _i respectively, which will provide us with some informa- 
tion about certain invariants of the system ([I]), namely, the differentiation index and the 
regularity of the Hilbert-Kolchin function. 

We introduce a double sequence [i/.^ of non-negative integers associated with the ma- 
trices Zk,i- 

Definition 5 For & € No and i € N> e _i, we define /j,k,i & No as follows: 
- fi Q ,i ■= 0, for every i € N> e _i; 
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- fjb kii := dimKker(3^) = k(n + r) - rank K (3fc,i), for k > 1 and i G N . 



Now, we will focus on the study of certain stationarity properties of the sequence 
((J>k,i)k,i- We begin by analyzing the behavior of the sequence (//fe,i)fc when the index i is 
fixed, which will be done by comparing the matrices Zk,i f° r & G N. 

First, let us observe the following obvious recursive relation which holds for every 
k > 1: 




Qjp(i-e+k+l) Qp(i-e+k+l) 
QQp-e+k + 1) QQ(i-e+k + l) 



dX( i + 1 '> 



QX( i + k + 1 ) 










gp(i-e+k + l) Qp{i-e+k+l) 



gj7(i+fc+l) 
QQ{i-e + k + l) 

QU(i+k+l) 



\ 



(2) 



When the differential system (TjQ) is linear (for instance, if k := Q and the system (JT]) 
is of type .At/ + BU = Y, where A and B belong to Q rxm ), the matrices Zki have a nice 
Hankel-block type form, but this is not exactly our situation. However, there is a main 
relation arising from their underlying differential structure which enables us to study them 
also in our non linear setting: 

Proposition 6 Let Z := Z±, . . . ,Z a be differential independent variables and let H be a 
differential polynomial in k{Z}. For all l,j G No the following relation holds in K: 

( 8H® Y dH^ dH® 

(3) 



\dZ^ +1 ) J dZd +1 ) dZW 

In particular, if H G {F, G} and Z G {X, U}, we have that Q£{j+\) = ^ ^ or ever y 3 — ^ + e ; 
since the order of H^ 1 ' is at most I + e, and therefore identity implies that 

g H (l+i) dH (i) 

Vj>/ + e. (4) 



dZd+V dZU) 

Note that these identities are also valid over the differential field IK, where the derivation 
is now the one induced by Xj = fj (X, U) in K. and all partial derivatives are regarded as 
elements in K. 

Due to the triangular form of Zk,i an d condition Q), all the matrices Zk,i have the 
same {n + r) x (n + m) — block in their main diagonal. 

Proof. Straightforward from the definitions and the Chain Rule. ■ 

We are now ready to prove the first stationarity property of the sequence {^k,i)k,i ( see 
also Remark H] above) : 
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Proposition 7 For each fixed i G N> e _i, the sequence (nk,i)keN ^ s non- decreasing and 
verifies the inequality 

r 

nmin{&;, e — 1} + minjfc, e — ej} < fik,i < min{A;, e} (n + r). (5) 
j=i 

In particular, there exists k G No (depending on i), 0</c<e + ra + Sj=i e j> suc/t i/iai 

Proof. The fact that (fJ*ki)k is a non-decreasing sequence follows immediately by observing 
that ker(3^ J x {0} C ker(^ +1 J for every fe G N (see ©, also Remark HJ. 

For every non-negative integer k G No, the matrix Zk,i has A;(?i + r) rows. Therefore, 
dimker(5^, J < k{n + r). On the other hand, due to Proposition [21 Remark [3] and the 
definition of /ik,i, we have that trdeg L (Frac( J 4i/Aj_ e+ i + fen J 4i)) = (m— r)(i + l)+e(n+r) — 
fi kji . Now, trdeg L (Frac(A i /A i _ e+1+fc fl^)) > trdeg L (Prac(^4 i /A D^)), since Aj_ e+1+fe n 
Ai C A n Ai, and so, the fact that the differential dimension of A is m — r implies that 
trdeg L (Frac(Aj/A n Ai)) > (m — r)(i + 1). Hence, < e(n + r) holds. 

In order to show the other inequality, we observe that, since the order of the polynomial 
gJ p) is p + ej (1 < j < r), the partial derivatives dG^ /dX® and dG^/dU^ are all 
zero for q > p + ej. In particular, we conclude that for t, s G N, with t > s, the derivatives 
<9G^ e+t) ^ ,^-(t+s) an j e+t ^ /Qjj( t + S ) are zero if i + s > i — e + 1 + ej, or equivalently, 
if e — ej > t — s. So, each polynomial Gj induces min{/c,e — ej} many null rows in the 
matrix Zj-i- Analogously, each polynomial Fj induces min{£;, e — 1} many null rows in this 
matrix. 

We conclude that the matrix Zk,i has at least nmin{£;,e — 1} + ^^ =1 min{/c,e — ej} 
null rows. Thus, the dimension of the kernel of the transpose matrix Z\ j (i-e. fik i) is at 
least nmin{A:,e — 1} + X^^=i rnin{fc, e — ej}. 

The second assertion follows directly from the fact that for every k > e, the inequality 
© reads n(e - 1) + Sj=i( e ~ e j) — A*fc,i < e ( n + r )- ■ 

In fact, in Theorem[9]we are able to prove a more precise result than that of Proposition 
[2 the sequence (/J>k,i)k€N is strictly increasing up to a certain index k{ < e + n + Y7j=i e j 
where it becomes stationary. 

For the sake of simplicity, we will use the following notations: 

Notation 8 The variables X, U involved in system fip are renamed in the following way: 
Zj := Xj for j = 1, . . . , n, and Z n+ j := Uj for j = 1, . . . , m (and the same is done for their 
corresponding formal derivatives). Analogously, the polynomials are renamed as: Hj := Fj 
for j = 1, . . . , n, and H n+j := Gj for j = 1, . . . , r. 

With these notations, the matrix 2k,i involves exactly the derivatives of the polynomials 
.ffO-e+i+p) with respect to the variables Z^- l+q \ with p = 0, . . . , k — 1 and q = 1, . . . , k. 

Theorem 9 Fix an index i G N> e _i and let ki G No be the minimum of all the k 's in 
No such that Hk+i,i = Hk,i (this minimum is well defined due to Proposition^. Then 
fJ-k,i = Hki,i for every k>k { . 



10 



Proof. According to Notation [8] we will rename variables and equations as Z := (X,U) 
and H := (F, G). 

The result is clear for ki = 0: in this case, fin = 0, which is equivalent to the fact 

that the matrix = qz(i+i) ^ as row ran ^' From identity (JH) in Proposition [6] 
we conclude that Zk,i has full row rank too or, equivalently, that fik,i = for all k. 

Now, let us assume that ki > 1. It suffices to show that the equality j = fJ-k-i,i f° r 
an arbitrary index k > 2, implies fik+n = f^ki- 

In the sequel, for a vector v E K i ( n + r ) we will write its description as a block vector 
v = (v±, ... ,v{) with Vj € K n+r . 

Due to the recursive relation ([2j), the identity ker(^ { ) x {0} = ker(^ +1 Jnjffc+i = 0} 
holds in K ( - k+1 ^ n+r ^> for every k € N and so, the equality fik,i = fi>k+i,i is equivalent to the 
inclusion ker(3^ +1 j) C {vk+i = 0}. Then, the theorem is a consequence of the following 
recursive principle: 

Claim: For all k € N, ker(^ J C {v^ = 0} implies ker(^ +1 J C {ffc+i = 0}. 

Proof of the Claim.- Let us show that if (vi, . . . , Ufe+i) € ker(^ +li ) then, the vector 

if = (wi, . . . ,Wk) € IK' c ( n+r ') defined as 



7***5 5 



tu fc := •Uj.+l, Wj := v j+ i - ibj+i, j = k-l 
lies in ker(^fi •), which implies the Claim. 

Our assumptions on the order of the equations imply that = ^ ^ or ^ ^ 

Hence, we have that w E ker(3^ A if and only if the following identities hold over K for 
every 1 < j < k: 

* Q H (i-e+t) 

We will proceed recursively for j = k, k — 1, . . . , 1. For j = k, the definition of w; and 
identity imply that 

QJj(i-e+k) Qjj(i-e+k+l) 

Wk dz(i+ k ) = Vk+1 dzt+ k +v = °' 

* Q H (i-e+t) 

Now, assume that > — — , . = 0. Differentiating this identity in K. and using 

t=j+i % 3 

identity (|3|) we get: 

A dH^~ e+ ^ -A / dH (i-e+£+l) dH (i- e +l)\ 

2- ^ SZ(*+i+i) + ^ 9Z(^+i) dZ^JT = °- 

£=j+i *?=j+i v / 
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This implies that 

dH^~ e+t > «9if(*- e +j) JL, Qff(i-e+£) * Qjj{i-e+l+l) 

i=j £=j+i e=j+i 

= ^(m + dz(i+j+1) +w k dz{ . +j+1) 

£=j+l 

k+1 Q H (i-e+t) Q H {i-e+l.) 
l=j+l 1=1 

where the second equality follows from identity (jH) and the third one is simply the definition 
of w. This concludes the proof of the theorem. ■ 

So far, we have studied the behavior of the sequence (fJ,k,i)k for an arbitrary (but fixed) 
index i £ N> e _i. In the remaining part of the section we will analyze the sequence (nk,i)i 
fixing the index k 6 N. 

We start by exhibiting a (non EC-linear) bijection between the kernels of the matrices 
$k,i and 2k,i+i for any index k £ N. 

Lemma 10 Let (v\, . . . ,Vk) and (w±, . . . ,Wk) be arbitrary elements in K fc ( n+m ) (here Vj 
and Wj denote vectors in K n+m for every index j). Then, for each i € N> e _i, the function 
9 : ]K fe ( n + m ) — > K fc(n+m ) defined as 9(v\, . . . , Vk) = v 2 - vi,v 3 -v 2 , ■ ■ ■ , Ufc - Vk-i) maps 
ker(3fc,?.+i) into ker(3fc,i). Moreover, 9 is a bijection between ker(3fc,j+i) and kev(Zk,i), mi/i 
inverse 9~ 1 (wi, . . . ,w k ) = (wi,w 2 + u>i,iu 3 + w 2 + u;^, . . . ,w k + w k -i H h tof* )■ 

Proof. It is easy to see that 9 is a bijection in ~K k ( n+m ^ with the inverse given in the 
statement of the Lemma. Let us show that it maps ker(3fc,i+i) to ker(3fc,i). 
We keep the notations introduced in Notation [HJ 

For arbitrary vectors (v%, . . . , Vk) and (w\, ... , tOfc) in 'K k ( n+m \ consider the following 
two families of sums (p = 0, . . . , k — 1): 



Ei+i v (v) := > — - — ., Vj and Ei v (w) := > 



it;,-. 



Note that v £ ker(3fc,i+i) if and only if Ei + i tP (v) = for p = 0, . . . , k — 1 and u> £ ker(3fc,i) 
if and only if Ei jP (w) = for p = 0, . . . , k — 1. 

First, we compare the vectors Zk,i+i v and Zk,i @( v ) fo r a given vector v = (v\, . . . , v^) G 

Qfj(i-e+2) QJj{i-e+l) 

K k(n+m) let observe that (f ) = £j o(f ) = #i o(#M), since - = - ,. , 

Qjj(i-e+2) Qff(i-e+l) 

for every j > 2 (recall that the order of H is smaller than e), .. „, = ,. 

(from ([3])) and 9{v)\ = v±. Now, for p > 0, we have: 
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Qff(i-e+2+p) f Qff{i~e+l+p) \ "\ ^ r^(i- e +l+p) 




Q Z (i+l+j) I Q Z (i+l+j) \ \ 3 ^(i+l+j) 



3=1 3=2 
= dZ«+l) "l+L dZ (i+3) fa-Vl) = ^W), 

where the second equality follows from identity ([3]), and the third one from the fact that 

— — — ,. ■ , — t = for p < k — 1. 

These equalities imply straightforwardly that 8 maps ker(3fc,i+i) into ker(3fc,i). 

In order to prove that it is onto, we may argue recursively: if w £ ker(3fc,i), then 
Ei >p (w) = for p = 0, . . . , k — 1. Now, £ , j + i j o(6 , ~ 1 (w)) = £'j + i i o(w) = Ei t o(w) = 0. Assum- 
ing that £'i_|_i i p_i(0 _1 (w)) = has already been proved, we deduce that Ei + i jP (8~ 1 (w)) = 
Ei <p (w) + (Ei + i^i(0~ 1 (w))' also equals 0. We conclude that 8~ 1 (w) € ker(3fc,i+i). ■ 

Even though the bijection 8 between ker(3fc 5 i+i) and ker(3fc,i) shown in the previous 
Lemma is not a K-linear map, it enables us to prove the following: 

Proposition 11 Let A; € No and i £ N> e _i be arbitrary integers. Then fiki = fik,i+i- 

Proof. The result is immediate if k = 0. Now, let k £ N be a positive integer. In order 
to prove that dimj^kerQ^ J) = dimic(ker(^ it suffices to show that and 3fc,i+i 

have the same rank, since they are two matrices of the same size. 
For each pair of indices j,t, 1 < j < k and 1 < t < n + m, set 

QJj[i— e+l,i— e+k] Qj^[i-e+2,i~e+k+l] 
Cjt '■= T^r^ and D»t : = j— — ., 

J ' dz[ l+l) h dz[ l+l+]) 

for the corresponding columns of the matrices Zk,i and 2k,i+ii respectively. 

Assume that a column Dj t of the matrix Zki+l i s a K-linear combination of the 
columns Djj to its right. Then, there exist elements a,j t S IK such that 

n+m k n+m 

D 30,to = Yl a 30,t D 30,t + Yl Y a 3,t D 3,ti 

t=t +i j=io+i *=i 

and so, the vector v := ( , . . . , , — a J0 , — ay +i> • • • > — a k) £ K fc ( n+m ) belongs to the ker- 
nel of Zk,i+i-> where denotes the null vector in K n+m , atj := (0, . . . , 0, 1, ctj 0t t +i, . . . , a JOi fe) 
and ctj := (a^i, . . . ,ctj t k) for j > jo + 1. By Lemma fTOl the vector 8(v) belongs to the 
kernel of Zki and, due to the particular form of the application 8, it turns out that the 
column Cj j is a K-linear combination of the columns to its right. Hence, the rank of 
is lower than or equal to the rank of 3fc,i+l- 

By means of the inverse map one proves in an analogous way that the rank of 
Zki+x is lower than or equal to the rank of Zki ■ ■ 
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Remark 12 The following alternative proof of this proposition was kindly suggested to us 
by Prof. F. Ollivier. From \14\ Prop. 10, Ch. IV] the dimensions of the kernels of2k,i an d 
3fe,i+i are the differential dimensions of the ideals generated in a differential polynomial 
ring in k{n + m) new variables by the linear equations defined by these matrices. Now, 
Proposition [77] follows since the bisection 9 induces a differential isomorphism between 
these fields. 

The previous proposition states that the sequence fik,i does not depend on the index 
i; therefore: 

Notation 13 In the sequel, we will write fi^ instead of n^ ii f or an U * £ N>e-i- 

So, Theorem [9] can be restated as follows: 

Corollary 14 There exists a non-negative integer a such that fik < Hk+i for every k < a 
and Hk = Mfe+i for every k > a. ■ 

4 The differentiation index 

Here we apply the properties of the matrices Zki established in the previous section to 
the study of well-known invariants of our DAE system (pQ): the differentiation index, the 
regularity of the Hilbert-Kolchin function and the order of its associated differential ideal. 

We start by introducing the notion of differentiation index of the DAE system (PQ) (see 
also |2], [8], [3], [20]): 

Definition 15 The non negative integer a introduced in Corollary^T^is called the differ- 
entiation index of the system |7p. 

Corollary 1141 states that the differentiation index a is the smallest non-negative integer 
k where the sequence /ij. becomes stationary. Since [i a = Hk for all k > a, we deduce from 
Proposition [71 

Remark 16 The following inequalities hold: 

r 

(e - l)n + y~](e - ej) < fi a < e(n + r); (6) 
i=i 

< o < min < e(n + r) , e + n + > . (7) 

For the last inequality see also [301 Section 5.2]. 
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4.1 The manifold of constraints 

For every i E N> e _i, the differentiation index a is strongly related to the minimum number 
of derivatives of the system (pTJ) required to obtain the intersection of the whole differential 
ideal A with the polynomial ring Ai, namely those polynomials in A which involve only 
derivatives up to order i (a similar result can also be obtained by rewriting techniques by 
means of [3Q|, Theorem 27]; see also [3 Lemma 14] and |21} Lemma 9]): 

Theorem 17 Let a E No he the differentiation index of the system {TJj. Then, for every 
i 6 N> e _i, the equality of ideals 

Ai_ e+ i +(T n Ai = A n Ai 

holds in the polynomial ring Ai. In particular, from inequality {?]) it follows that the 
identity A i+1+n+ Y^ r _ i ej D Ai = A D Ai holds for every i 6 N> e _i. Furthermore, for every 
i € N> e -i, the differentiation index a verifies: a = min{/i € No : Aj_ e+ i + ^n^4i = AnAi}. 
Taking i = e — 1, we have 

cr = min{/i G M : Afc fl i e _i = A n i e _i}. (8) 

Proof. Fix an index i € N> e _i. Let us consider the increasing chain (Aj_ e+ i + fc n^4j)fc g pj 
of prime ideals in the polynomial ring A^. From Proposition Remark[3]and the definition 
of the sequence [i^, for every k € No, we have 

trdeg L (Frac(A i /A i _ e+1+fc n A)) = (ro - r)(i + 1) + e(n + r) - /i fe . (9) 

Since fi^ is stationary for k > a (Corollary |14|). the previous equality implies that all the 
prime ideals Aj_ e+ i + fc n Ai have the same dimension for k > a, and so, they coincide 
because they form an increasing chain of ideals. On the other hand, any finite system of 
generators of the prime ideal A n Ai C Ai belongs to Aj_ e+ i + fc n Ai for all k big enough, 
which finishes the proof of the first assertion of the Theorem. 

In order to prove the second part of the statement, for each i £ N> e _i, let hi be the 
smallest non- negative integer such that Aj_ e+1+ / li n Ai = A n A{. By the definition of hi, 
the transcendence degrees trdeg L (Frac(^4j/Aj_ e+ i + fc n Ai)) coincide for k > hi, and so, 
is constant for k > hi (see identity ((HJ) above). This implies that a < hi. The equality 
follows from the first part of the statement and the minimality of hi . ■ 

Identity (jHJ) in Theorem [T7| can be regarded as an alternative definition of the differ- 
entiation index (see also |21^ Section 3.2] for first-order DAE systems). In particular, it 
gives the following interpretation of the differentiation index a (see [8] for the case e = 1): 

Remark 18 If the differentiation index a of system (CP equals 0, there are no constraints 
on initial conditions for the system (recall that Aq := (0)). In the case when a > 1, the 
quantity a — 1 is the minimal number of derivatives of the equations in the system needed 
to obtain all the relations that the initial conditions must satisfy (the so called "manifold 
of constraints on initial conditions"). 
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Another fundamental property of the differentiation index, concerning the number of 
derivatives of the equations required to obtain an explicit equivalent ODE system, will be 
considered in Subsection 14.2.31 and Section [5] below. 

First we will show how our previous results can be applied in order to estimate the 
regularity of the Hilbert-Kolchin function. 

4.2 The regularity of the Hilbert-Kolchin function and applications 

We recall the main basic facts of the Hilbert-Kolchin function for the particular case of 
our input system ([I]) and our ideal A (see [141 Chapter II] for the general theory). 

The Hilbert-Kolchin function : No — ► No of the differential ideal A C L{X, U} is 
defined as 

H A (i) := trdeg L (Frac(A i /(A n A { ))) 

for every i € No- Since the ideal A we are considering has differential dimension m — r, 
the identity H&(i) = (m — r)(i + 1) + ordi,(A) holds for i sufficiently big (see for instance 
[141 Ch. II, Sec. 12, Th. 6]), where ordiL,(A) is a non-negative integer depending only on 
the differential ideal A which is called the order of the ideal. The polynomial 

H A (T) := (m - r)(T + 1) + ord L (A) 

is called the Hilbert-Kolchin polynomial associated to the ideal A. The regularity of the 
Hilbert-Kolchin function H/± is defined to be the minimum integer iq £ No such that 
Ha{}) = T^aW holds for all i > %q. 

4.2.1 The regularity of the Hilbert-Kolchin function 

The theory of characteristic sets (see for instance [28]) can be used to give a precise 
estimation of the regularity of Ha- let C be a characteristic set of the differential ideal 
A C L{X, U} for an orderly ranking (see [28] ) in the variables X, U. Then, the regularity 
of H A is equal to max{ord(C) : C G C} — 1 (this fact follows as an immediate consequence 
of the proof of p3J Ch. II, Section 12, Th. 6 (d)], see also [1 Theorem 3.3]). 

The results developed so far enable us to exhibit the following simple upper bound for 
this regularity, which depends only on the order of the polynomials involved in the system 

©: 

Theorem 19 The regularity of the Hilbert-Kolchin function of the ideal A over the ground 
field L is bounded by e — 1. 

In particular, for first- order systems of type {7]) (in other words, for the case e = 1), the 
Hilbert-Kolchin function of A coincides with the associated Hilbert-Kolchin polynomial. 

Proof. It suffices to show that the equality H^(i + 1) = Ha(i) + (m — r) holds for every 
i G N> e _i. 

Fix an index i € N> e _i. Due to Theorem 1171 we have that A n Ai = Aj + i_ e+cr n A,- t 
and A n A i+1 = A i+2 - e +a n A i+1 . So, H A (i + 1) = trdeg L ( J 4i + i/(A i+2 _ e+CT n A i+1 )) and 
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H&(i) = trdeg L ( J 4j/(Aj + i_ e+CT C\Ai)). Thus, using Proposition [2] (recall also Notation [T3|) 
we obtain: 

H A (i + l) = (m-r)(i + 2) + e(n + r)-fi a , 
H&(i) = (m — r)(i + 1) + e(n + r) — 

Hence, the equality H&(i + 1) = H&(i) + (m — r) holds. ■ 

As we have already pointed out, the regularity of the Hilbert-Kolchin function can 
be obtained from the maximal order appearing in any characteristic set associated to 
the differential system (TjQ) for an orderly ranking. Therefore, our upper bound on the 
regularity implies an upper bound for this maximum: 

Corollary 20 Let C be a characteristic set of A C L{ V, U } with respect to an orderly 
ranking on the variables X, U . Then, the order of each of the differential polynomials in 
C is bounded by the maximal order e := max{l, e±,e2, ■ ■ ■ ,e r }. ■ 

4.2.2 Upper bounds for the order of the ideal A 

The proof of Theorem [TU] also shows an indirect relation between the Hilbert-Kolchin 
polynomial and the differentiation index a. More precisely: 

Remark 21 The Hilbert-Kolchin polynomial Ti&.(T) of the ideal A can be written as 
T~(-a{T) = (m — r)(T + 1) + e(n + r) — fx a . Equivalently, ordL(A) = e(n + r) — fj, a . 

From inequality ([6]) and Remark [21] we deduce the following upper bound for the order 
of the ideal A in terms of the orders of the defining equations: 

Corollary 22 The order of the ideal A can be bounded as follows: 

r 

ord L (A) < n + ^ej. ■ 
j"=i 

This result is known as Greenspan's bound. It was first proved by B. Greenspan in [S] 
for difference equations and extended later by R. Cohn in [6] for O-dimensional components 
of differential systems. 

Another upper bound for the order is given by Ritt in [28L Ch. VII, p. 135] (see also [14^ 
Chapter IV, Proposition 9] and [31j ) and holds for general O-dimensional DAE systems: 
the order of its associated differential ideal is always bounded by the sum of the maxima 
of the orders of each variable in the given differential polynomials. This bound is called 
by Ritt a differential analogue to Bezout's Theorem. In our case, we are able to give a 
simple proof of the following extension of this fact: 

Corollary 23 For each i = 1, . . . m set := maxjord^Gj, j = 1, . . . , r}. We have the 
following inequality: 

m 

ord L (A) <n + y~]ej. 

i=l 
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Proof. We transform our input system ([I]) in an equivalent first-order system by introduc- 
ing new variables and equations in the usual way: for each i = 1, . . . , m and £ = 0, . . . , e» 

(£) 

we define Ua := U- . Making this substitution in the original system and adding the 
equations Ui t e = Ui,t—\ (1 < i < m, 1 < I < ei) we obtain a first-order DAE system 
(E') with a structure similar to the original system ([T]), but such that the order of each 
parametric equation gj = Yj is (in other words, the new e,-'s are zero for j = 1, . . . , r). 

Set U for the set of variables Ua with 1 < i < m, < £ < e{. For each integer p > 0, 
denote by A' C k{X,U} the differential ideal generated by the equations defining (E') and 
by A' p the polynomial ring h[X^,U^}. 

It is easy to see that for each p > e the map 

u(j) ^ f ^ if < j < 

induces a monomorphism A p / (A n A'/ (A' n vl' ). Hence, for p sufficiently big, the 

inequality Ha(p) < T~(-A'(p) holds. Since the differential dimension of both systems (pQ) 
and (E') is m - r, we conclude that ordL(A) < ordL(A'). 

The proof finishes by applying Corollary|22lto the first-order system (E'), and observing 
that in this case the bound is simply the number of "non parametric" equations (since the 
new parametric ones have order zero), which is exactly n + YliLi e «- H 

An interesting and more precise Bezout-type bound for 0-dimensional systems (which 
includes both Corollaries [22] and f23|) is discussed by Jacobi in [12] (see also [211 Ch.VII, 
§6]). Even though Jacobi's bound remains conjectural for the general case, in |15] (based 
on [29]) it has been proved to hold for a class of 0-dimensional differential ideals including 
those considered in this paper. More precisely, suppose m = r and set e,j for the maximum 
of the orders of a derivative of the variable £/$ occurring in gj, putting := — oo if the 
variable Ui does not appear in gj. Then, Jacobi's inequality states that 

m 

ordjL(A) < n + max e T (j)j, 

where S m is the permutation group of order m. Now we show an example where Jacobi's 
bound is reached and improves the previous ones. 

Example 24 Consider the following first-order DAE system with coefficients in k := Q: 

U 2 + Ui 
u m + U m -i 

with m > 2. Here, both upper bounds for ordt(A) given by Corollaries HH and\2^ equal 
m — 1, while Jacobi's upper bound is 0. It is not difficult to see that ordi > (A) = 0. 
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4.2.3 The index and an equivalent explicit ODE system 

The estimation for the regularity of the Hilbert-Kolchin function allows us to give a first 
result concerning the number of derivatives of the input equations required to obtain an 
equivalent explicit ODE system in the O-dimensional case. We will show that this number 
is at most the differentiation index of the system: 

Theorem 25 Let (£) be a DAE system as in |7P of differential dimension (or equiv- 
alently, r = m), maximal order bounded by e and differentiation index a. Let S = 
{£l , ...,£s } C {Xl 6 " 1 ! , C/"^ 1 !} be an algebraic transcendence basis of the fraction field 
of A e _i/ A n A e _i over the field L. Then: 

1. for each i = 1, . . . , s there exists a non-zero separable polynomial Pi with coefficients 
in the base field k, such that ?j(yW,S,^ e) ) G (pH^M) C k[Y^,X^+^,U^ +e i]; 

2. set {773+1, . . . ,r] n+m } := {X,U} \ {£1, . . .,£,}. Then, for all i = s + 1, . . . ,n + m, 
there exists a non-zero separable polynomial Pi with coefficients in the base field k, 
such that P^Y^-^^rf^) G {F^~ x \ G^~^) C k[Y^~ l \X^\ U^+ e ~\ 

In particular, for every 77 G {X, U} there exists a separable non-trivial polynomial 
relation between rfi e ' and E modulo A which can be obtained using at most a derivations 
of the input equations. 

Proof. In order to prove the first statement, we fix a variable £j . Since we are in a 
differential O-dimensional situation, from the upper bound on the regularity of the Hilbert- 
Kolchin function (Theorem 1 19j) . we have that the set E is also an algebraic transcendence 
basis of the fraction field of the ring A e / A n A e . Then, there exists a polynomial Q in s + 1 
variables with coefficients in L, such that Q(3, Q^) belongs to the ideal An^e = A a+ iPiA e 

° dQ 

(Theorem 1 1 Th . Clearly, this polynomial can be chosen so that — p-r 4. A (i.e. separable). 

de e > 

^0 

Now, the polynomial Pj of the statement can be easily obtained from the previous Q by 
multiplying it by an adequate factor in &{Y} and evaluating the superfluous variables Y^ 
for all I > a + 1 at suitably chosen elements of the base field k. 

The second assertion follows similarly, but in this case we use the fact that the family 
{77^ ^ , S} is algebraically dependent when regarded in the fraction field of -A e -i/ A n A e ^\ 
over L, for all % = s + 1, . . . , n + m. ■ 

The next section is devoted to showing slightly more precise results in the same spirit 
of the previous Theorem. 

5 Toward an equivalent explicit ODE system 

In this section we make use of Theorem in order to give more precise results on the 
number of derivations required to obtain from a system of type ([TJ an equivalent explicit 
ODE system. In the first subsection, we show that for O-dimensional systems there exists 
a more accurate upper bound for the number of derivatives and the order of the involved 
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variables by means of an adequate modification of the notion of differentiation index. In 
the second subsection, we consider the case of positive dimension for first-order systems. 
Finally, we discuss some quantitative and symbolic algorithmic aspects. 



5.1 Equivalent explicit ODE systems and an alternative notion of dif- 
ferentiation index for zero-dimensional systems 

Throughout this subsection, our input system defined in (TjQ) is considered to be O-dimensional 
and further, it is assumed that no variables X appear (i.e. r = m and n = 0). More ex- 
plicitly, we consider the following DAE system with generic second members: 

' 9l (U,U,...,U^) = *i 

(£) ■= \ : (10) 

k g r (U,U,...,U^) = Y r 

where U := U\, . . . , U r and Y := Y\, . . . ,Y r are sets of differential indeterminates, e,- S No 
for every 1 < j < r, and e := max{ei, . . . , e r }. We suppose e > 1. 

For each 1 < i < r we set 6j for the maximum of the orders of the variable £7j in the 
polynomials g\, . . . ,g r . Since we assume a differential O-dimensional situation, all variables 
must appear in the system; that is, we have that < q < e for all i. 

We introduce a new set of differential variables Z := Z±, . . . , Z r verifying: 

Ui = z\ e - ti) (11) 

for i = 1, . . . , r. Writing the input system (X) in terms of the indeterminates Z we obtain 
a new DAE system: 

' 9l (Z,Z,...,Z^) = Y x 
(E) := : (12) 

k g r (Z,Z,...,Z^) = Y r 

Note that each variable Z{ appears at order exactly e in at least one of the equation in 
the system (£), and e is the maximum of the orders of all variables appearing in (£). 

We denote a the differentiation index of the system (S). The number a can be regarded 
as an invariant of the DAE system (E) and it represents, in some sense, a more accurate 
version of the global differentiation index introduced above, allowing us to distinguish the 
variables more precisely according to their maximal orders (c.f. the notion of structural 
index introduced in [23] . [24] ) . We first show that the differentiation index a of (E) is 
always an upper bound for a: 



Lemma 26 The inequality a < a holds. 

Proof. Let G be the polynomials defining the system (E), and let A be the differential 
ideal generated by the differential polynomials G in the differential polynomial ring L{Z}. 
For each p € No, we set A p for the polynomial ring L[Z^] and A p for the ideal generated 
by the polynomials & p ~^ in the ring A p ^\ +e (note that the maximum of the orders of the 
polynomials G equals e). 
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Denote by Z the set of indeterminates Z\ where i = 1, . . . , r and < ii < e — ej (if 
e = fj no derivatives of appear in the set Z). The variables Z will be considered as 
algebraic independent variables with respect to the original variables U and their successive 
derivatives. We also use the following standard notation: if A is a ring, 21 C A an ideal 
and X a set of indeterminates over A, 2l[X] denotes the ideal of A[X] consisting of all the 
polynomials whose coefficients belong to 21. 

In order to prove the inequality a < <r, applying Theorem 1 171 to (£), it suffices to show 
that the inclusion of ideals 

A n 2 e _i c A CT n 2 e _i 

holds. 

Let H G A n A e _\. Since H G A, there exists a sufficiently big index q G N such that 
If G Aq. So, we can suppose that H is a polynomial combination involving only variables 

~{p) 

Z of order at most q of polynomials of type Gj with p < q. 

(e—e-\-i) (£) 

Replacing Z\ 1 with U- , for i = l,...,r and each non negative integer I, we 
see that the polynomial H lies in the ideal A 9 [Z] n L[C/J; ei_1 ', . . . , C/r er_1 ^][^] (note that 

L[?7j £1_1] , . . . , ul er ~ l] ] [Z] = 2 e _i and that if Ci = neither U { nor its successive derivatives 
appear in this polynomial ring). Since e% < e for all i, we conclude that H G A 9 [2] n 
A e -i[Z] = (A q n A e _\)[Z\. Now, due to identity (JSJ) in the statement of Theorem 1171 we 
have that A q n A e _i C A CT n A e _i. We infer that H G (A CT n A e _i)[2] = A CT [Z] n A e -i[Z}. 
The result follows observing that A CT [2] = A CT and j4 e _i[^] = A e _i. ■ 

The following pendulum-type system shows that the inequality in the previous Lemma 
may be strict: 

Example Consider the DAE system: 

Y x 

Yi ■ 

Y 3 

The corresponding matrix 5s,i (see Definition [T]) is 



[ U^ + UxUs = 
(E) := U^ + U 2 U 3 = 
I vf + v*-! = 
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3f7i 


U 3 
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3f/ 3 


3t> 2 





^3 


U-2 
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6J7i 


6U 2 





2E7i 


2C/ 2 

































617, 



(2) 



\ 12U 



(2) 



6U. 



(2) 



6U- 
6U. 



(2) 
L 

(2) 



12U, 



(2) 



4U 3 4C/i 
4Us 4t/ 2 
8Ui 8U 2 



Us C7i 
Us u 2 
2Ui 2U 2 









1 



1 
0/ 



The dimension (over the field k(U)) of the corresponding kernels of Z\ i are \ii = i, for 
: 0, ... ,4, and (1$ = 4. Therefore, the differentiation index of the system (E) is a = 4. 



Changing the variables as in (jTTj) : Z\ := U\, Z 2 
the previous system into: 



(£) :- 



Z 



(2) 
1 

,(2) 



(2) 



(2) 



Z 2 + z 2 z^ 
z\ + zl-\ 



U 2 and Z. 



Ys 



(2) 



Us, we transform 



Here, the corresponding Jacobian sub-matrix ^31 is: 



/ 1 













1 


z 2 























Zi 


1 








z 2 

















7(2) 





z? 








7(2) 

Z 3 


7(2) 





V 2Z X 


2Z 2 













Zi 

z 2 



















\ 










2Zi 1 Z\ 
2Z 2 1 Z 2 
/ 



and the dimensions of the kernels of the matrices Z\ \ (for k = 0, 1, 2, 3) are /j,q = 0, fj,\ = 1, 

Hi = 2 = fj,s- Therefore, the differentiation index of the system (E) is a = 2, which is 
smaller than the index a = 4 of the original system (E). 
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Remark 27 Applying a standard change of variables to the system (£) defined in U0\) 
we obtain an equivalent first- order system (£) with differentiation index a. With a similar 
argument to that of Lemma\2^ it is easy to see that a < a < a. In the previous example, 
we have that a = 3, which shows that the inequalities may be strict. 

By means of the modified differentiation index a, it is possible to give a more accurate 
version of Theorem [25} 

Theorem 28 Let (£) be a DAE system as in (Q7J). For i = 1, ...,r we denote ej the 
maximum of the orders of the variable Ui in the equations defining (£). Then there exists 
a set hi C {u[ ei } ...,Ur } (here, if €{ = 0, no variable U± appears) and separable 
polynomials Pi, . . . ,P r such that 

Pi(Y^,U,U^) G (Gf\...,G^). 

In other words, there exist non-trivial (separable) polynomial relations between each U-^ 

and a fixed family of derivatives U C {u[ ei ^, . . . , uY T } (which is a family of order 
strictly lower than maxjje.;} — 1 ) that can be obtained from the first a many derivatives of 
the input equations. 

Proof. Let (£) be the DAE system obtained from (£) after the change of variables (1111) . 
For j = 1, . . . , r, denote by Gj the polynomial obtained from Gj after this change. 

Applying Theorem [25] to (E), there exists a subset 5 C {z| e , . . . , zl e ^} such that 
for i = 1, . . . , r there exists a separable polynomial Qi satisfying 

Q t {Y^\~,z\ e) ) e C k[Y^,Z^ + % (13) 

Fix an index i and its polynomial Qi. Writing ujp for each variable zj e ej+ ^^ we have 

that the set S may be decomposed as S = Ejj U where Ejj C {lf[ ei , . . . , lll €r ^} 

and E z C {Z { *~ ei ~ l \ . . . , zl e ~ €r ~ 1] }. With this notation, condition (035 can be written 

as Qi(Y^\ Ejj, Ez, U- e ^) G (G^), where the ideal (G^) is considered in the polynomial 

ring k[YW,Z [ r ei ~ 1] , ■ ■ ■ , 4 e ~ £r-1] , Uf +ei \ . . . , uf +e \ 

Then, the Theorem follows from (|13p taking U := Ejj and evaluating the variables 

zj e ej for j = 1, . . . , r in generic elements of /c. ■ 

We remark that Theorem [28l improves upon Theorem 1251 in at least two points: first, 
the number of derivatives required to obtain the polynomial relations is the modified 
differentiation index a instead of a. On the other hand, the order of each variable Ui in 
the corresponding relation Pj is the maximum order of the single variable Ui in the system 
(namely, ej) instead of the maximum of the orders of all variables (namely, e). 

5.2 Differential transcendence bases and equivalent ODE form for generic 
first-order DAE systems of positive dimension 

In the previous sections (more precisely, in Theorems and [2"5|) we saw how to obtain 
implicit equations for derivatives of low order of each variable in a DAE system after at 
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most a (or a) many differentiations of the input equations. In both results a "square" 
O-dimensional situation was assumed. 

Here we show a result of the same kind in the case of positive differential dimension 
for generic first-order systems of type (pQ). Even though the straightforward idea could 
be applied (namely, the localization into a differential transcendence basis to reduce the 
problem to a O-dimensional situation), different transcendence bases may lead to systems 
with different quantitative properties. For instance, the number of derivatives required to 
obtain an equivalent explicit ODE system may change. In this sense, there exist transcen- 
dence bases which are more adequate for our global analysis. Before stating the precise 
results we consider a simple example in order to illustrate this kind of phenomena. 



5.2.1 Example 



Consider the following 1-dimensional first-order DAE system with coefficients in k :- 
borrowed essentially from Section 3.4]: 



y-2 



Y 



m—1 



Ui + U m 
U 2 + Ux 

Um-l + U m -2 



(14) 



with m > 2. 

Here, we have n 
3i,o G K( m - 1 ) xm is 



and r 



m 

( o 
1 



1. The associated field K is 



and the matrix 



1 \ 





1 0/ 



which has full row rank m — 1. Hence, we have (Mi = [Xq = and so, Definition 1151 states 
that the differentiation index a of the system (j!4j) is equal to 0. 

Following the slightly informal definition of the differentiation index given in the In- 
troduction, the fact that a = should imply that all the derivatives £/j can be written 
in terms of the variables U\, . . . , U m using the equations (i.e., the variables Y±, . . . , 7 m _i). 
While this is trivially true for all the variables different from U m —i, 



Ui 



Y x 
Y 2 



U 2 



(15) 



U„ 



Y 



m—1 



it is not so clear how to find a similar relation for U m -x, because it does not appear in the 
equations. 

It is quite natural to consider U m -\ as a free variable and then to interpret the previous 
system as a O-dimensional system over the field k\ := Q(U m -x}- Now, the matrix 3i o is 
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the (m — l)-square matrix 

/ ••• 1 \ 
1 ••• 

\0 ■•• 1 0/ 

which is clearly non singular. Hence jjL\ = jjlq = and so, the new associated differentiation 
index o\ is equal to 0. The relations (|15[) can be seen now as a full rewriting of the 
derivatives in terms of the variables (c.f. Theorem 125 j) . Furthermore, we observe that the 
order of the ideal associated to the system does not change by the extension of the ground 
field from Q to Q{U m -i): following Remark [2T| we have ord^u m _ 1 ^(A) = ordL(A) = m—1, 
since the sequence (fJ<k)k is the same in both cases. 

On the other hand, if in the system (|14p we take U m as the free variable and consider the 
original system as a (m — l)-square 0-dimensional DAE system over the ground differential 
field k 2 := Q{U m ), the matrix 3i,o is the (m — l)-square matrix 

/ \ 

1 ••• 

\0 ■■■ 1 0/ 

So, fj,i 7^ = no and therefore, the differentiation index o 2 of the system (|14jl over k 2 is 
strictly positive. In fact, it is easy to see that a 2 = m — 1. Now, Theorem 1251 ensures that 
the order of derivatives of the equations required to express the derivatives Ui, . . . , U m -x i n 
terms of U\, . . . , U m -\ is at most m — 1. Indeed, the variables U\, U2, • • • , U m -\ can be writ- 
ten in terms of derivatives of the equations (i.e. derivatives of the variables Y%, . . . , Y m _i) 
and elements of the base field hi as follows: 

U 2 = Y 2 -Y l + U ( ^> 

U 3 = Y 3 -Y 2 + Y{ 2) -U^ 

= y m -i-y m - 2 + --- + (-ir- 2 n (m ~ 2) + (-ir _1 ^ m " 1) - 

Observe also that the order of the ideal is not preserved after localization in U m , since 
ord L(l/m) (A) = f m - 1 = ord L (A). 

The previous example shows two different situations arising when considering different 
transcendence bases. The first localization of the system (|14p seems to follow the behavior 
of the original system more closely than the second one, since in that case the derivatives 
of the unknowns which are not in the transcendence bases could be written in terms of 
the unknowns themselves using as many derivatives of the equations as the differentiation 
index of the original system. 

In the next two subsections we will show that for any first-order system of type ([T]) there 
are suitable differential transcendence bases which enable us to obtain relations between 
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the remaining differential dependent variables using "few" (as many as the differentiation 
index) derivatives of the equations. 

5.2.2 Differential transcendence basis preserving the order 

In the sequel, we suppose that the input system (E) is of first order (or equivalently, 
e = 1). The notations correspond to those introduced in the previous sections. 

We will denote T% := Frac(^4j/(A n Ai)) for every i G No- As before, we will use the 
same notation for an element of Ai or its class in T% whenever the ring in which it is 
considered is clear from the context. The fact that T Ti+i for every i G No allows us 
to consider any subset of T as a subset of Tj for every j > i, which will also be done 
without changing notations. 

Lemma 29 Let B C Ai and let ( G Ai be a polynomial such that its class ( G T is 
algebraic ouerL(S). Then, £ G T+i is algebraic over L(S U S) , where B denotes the set 
of classes of all derivatives of elements in B. In particular, if L(S) T is an algebraic 
field extension, then L(S U B) <— > T+i is also algebraic. 

Proof. The result is immediate if £ G B. So, let us consider the case when Q ^ B. Let 
P G L(S)[T] be the minimal polynomial of £ with respect to the field extension L(S) Ti. 
Multiplying it by a non-zero element in L[S], we may assume that P G L[S,T] and has 
non-zero leading coefficient. 

We have P{B, () G A n A u and so P(B UB,(,() £ An A i+l . Now, P(B U S, (, () = 
Q(BUB, C) + §p C f° r some polynomial Q. As deg T (|y) < deg T (P), the minimality of 
P implies that §£(S, () $ A and so, P(SUS, (, T) is a non-zero polynomial in L(SuS, [T] 
annihilating £ in T+i. This implies that Q is algebraic over L(S UB, (). 

Since the field sub-extension L(S U B) L(S U B, () of L(S U B) T + ± is algebraic, 
we conclude that Q is algebraic over L(S U S). ■ 

Proposition 30 Let s := ordiL(A). There exists disjoint subsets W := {Wi, . . . ,W m - r } 
and£ := of the set {Xj, . . . , X n , Ui, . . . , U m } such thatBi := {W{ 1 \ . . . , W^L r , £i, 

is a transcendence basis of the algebraic field extension L "-^ T% for every i G No- In par- 
ticular, the set {Wi, . . . , W m - r } is a differential transcendence basis of the differential field 
extension L Frac(L{X, U}/A). 

Proof. Let So C {X±, . . . , X n , U\, . . . , U m } be a transcendence basis of L To- Then, 
L(So) > is an algebraic field extension and so, due to Lemma 1291 L(So U So) > T\ 
is also algebraic. Hence, Bq U So contains a transcendence basis of T\ over L. Since So is 
algebraically independent over L, and trdeg L (.Fi) = m — r + trdeg L (.Fo) (see the proof of 
Theorem [T9|) . there exists a subset So C Bq with m — r elements such that B\ := So U Bo 
is a transcendence basis of the extension L T\. 

Let us denote Wi , . . . , W m - r the variables whose first derivatives are all the elements 
in S (note that {Wi, . . . , W m - r } C S ) and let {6, •••,6} := B \ {Wi, . . . , W m _ r } 
(observe that #So = (m — r) + s since the Hilbert-Kolchin function coincide with its 
associated polynomial as shown in Theorem I19p . We will show that, for every i£N, the 

set Bi := {Wi, ■ ■ ■ , wJ^L r , £i, . . . , £ s } is a transcendence basis of L <-^> T. 
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The case when i = 1 follows from our previous construction. Let us assume now that Bi 
is a transcendence basis of L T% for a fixed positive integer i 6 N. Then, by Lemma [29l 
L(i3.;U£>j) = L(£>j + iU{£i, . . . , £ s }) ■— > is an algebraic field extension. Now, L(£>j + i) 
L(i3j + i U {£i,...,£ s }) is an algebraic sub-extension of L(i?j+i) ^i+i, since each of 
the elements £j is obviously algebraic as an element of L(Sj + i) over L(£>i). Therefore, 
L(Bj + i) is an algebraic extension and, taking into account that trdeg L (J-"j+i) 

equals the cardinality of Bj+i (see Theorem 1 19h . we conclude that Bj+i is a transcendence 
basis of L ■— > ^i+i- ■ 

It is clear from the proof that the differential transcendence basis {W±, . . . ,W m - r } 
preserves the order of the ideal (the constant term of the Hilbert-Kolchin polynomial) 
after extending the ground field L to L(Wi, . . . , W m - r ). 

5.2.3 An equivalent ODE system 

Theorem 1 191 and Proposition [3U] enable us to deduce the following "implicit function type" 
result in terms of the differentiation index a introduced in Definition [15] (see also [8] Section 
3]), an analogue of Theorem [25] for first-order non-square generic DAE systems: 

Corollary 31 The variables X, U can be split into three subsets W := {W±, . . . , W m - r }, 
£ := {£i, ... ,6} and rj := {rj s+ i, . . . ,rj n+r } so that: 

1. W is a differential transcendence basis o/L Frac(L{X, U}/A); 

2. U £ is an algebraic transcendence basis o/L J-j for all j £ No; 

3. for each i = 1, . . . , s there exists a non-zero separable polynomial Pi with coefficients 
in the base fieldk, such that Pi(Y^,W,W,^,ii) & (F^,G^) C k\Y^ ] X^ +1] U^+% 

4- for each i = s + 1, . . . , n + r, there exists a non-zero separable polynomial Pi with 
coefficients in the base field k, such that P^Y^- 1 } ,W,£,r)i) G (irk- 1 ), G^- 1 )) C 

Proof. Let W and £ be subsets of variables as in Proposition [30] Then, the first and 
second conditions in the statement hold. Let rj := {r] s+ i, . . . ,r] n+r } be the set of the 
remaining unknowns X, U (i.e., those different from the W 7 s and the £'s). 

Now the proof runs mutatis mutandis as in Theorem [25] For instance, for the fourth 
item we observe that for every i, s + l<i<n + r, the element rji £ J-q is algebraic over 
L(W U £), since W U £ is a transcendence basis of L ^ Tq. Hence, there is a non-zero 
polynomial Pi with coefficients in L such that Pi(W, £, rji) € An^4o = A a nAo ( see Theorem 
fT7|) . In the same way, for the third item we have that W U W U £ is a transcendence basis 
of L «Fi which ensures the existence of non-zero polynomials Pj (z = l,...,s) with 
coefficients in L verifying Pi(W, W, £, ^) G Anii = A a+1 n Ai. 

The polynomials Pj's that we are looking for can be easily obtained from the previous 
Pj's, by multiplying them by adequate factors in k{Y} and evaluating the superfluous 
variables Y^ at suitably chosen elements of the base field k (for all I > a in the case of 
i > s + 1 and for all / > a + 1 for i < s). ■ 
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Formally, the fourth statement of the previous Corollary makes no sense if the differen- 
tiation index a is zero. However, it admits a natural interpretation also in this case: this 
situation corresponds exactly to the case where the Theorem of Implicit Functions can be 
applied in order to write each derivative of the variables {X, U} \ {W} in terms of the 
same variables and derivatives of the VF's (obviously not in a polynomial nor a rational 
way). So, the fourth item must be empty. More precisely: 

Remark 32 Under the conditions of Corollary \31\ suppose also that the differentiation 
index a is zero. Following Remark \21l the Hilbert-Kolchin polynomial of A over L is 
%a(T) = (m — r)(T + 1) + (n + r) and the order of this ideal is n + r (recall that we 
assume e = 1 and fiQ is defined to be 0). So, the set n is empty, or equivalently £ = 
{X, U} \ {W}. Then, there exist non-zero polynomials Pi, i = 1, ...,n + r, such that 
Pi(Y, W, W, C ii) G (F, G) C k[Y, X, X, U, U}. 

5.3 Quantitative and algorithmic aspects 
5.3.1 Degree bounds of the implicit equations 

The non-zero polynomials Pi of Theorems [25] & [28] and Corollary [3T| are not uniquely 
determined without additional requirements (as minimality of order and degree, irre- 
ducibility, etc.). However, the conditions stated in those results allow us to choose a 
family of such polynomials that can be regarded as eliminating polynomials of suitable 
algebraic-geometric situations, which enables us to estimate their degrees. 

In order to illustrate these facts consider for instance the situation of Corollary [3T1 let 
A; be a fixed algebraic closure of the ground field k. For each N £ N we denote by A N 
the affine space k N equipped with the Zariski topology. Set Nq := (r + n + m)a + n + m 
and TVi := (r + n + m)(o~ + 1) + n + m and let Vo C A^ and Vi C A^ 1 be the algebraic 
varieties defined by the ideals (F^ , C fc^UM, J/M] and (F^,G^) C 

k[Y^\X^ +l \U^ +l ^\ respectively, that is: 

V := {F^- 1 ! = , G [a ' 1] = 0} and Vi := {F^ = , G [a] = 0}. (16) 

Note that both varieties are irreducible complete intersection and their dimensions are 
m(o~ + 1) + n and m(a + 2) + n respectively. 

Let W := {Wi, . . . ,W m - r }, £ := {£i, ... ,6} and n := {r] s+1 , . . . ,r] n+r } be a partition 
of the set of variables {X, U} as in Corollary EH We define linear projections 0i (i = 
1, . . . , s) and 7Tj {i = s + 1, . . . , n + r) as follows: 

0i : Vi -> A r ( fT_1)+2m+s+1 , 6i{y^\x^ +l \u^ + ^) := (yH ,w,w, 
-Ki : V -> A ra+m - r+s+1 , Ki{y^- l \x^\uW) := (y^~^ , w, £, Vi)- 

From Proposition [30l we deduce that the set {Y M , W, W, £} (resp. {Y^~^,W^}) 
is algebraically independent in the fraction field k(Vx) (resp. k(Yo)) over k. On the 
other hand, due to Corollary EH the set {yW, W, W, £, ii} (resp. {Y^~^ , W, £, ty}) is 
algebraically dependent. Thus, the closure of the image of the map 0i (resp. 7Tj) is a fc- 
definable irreducible hypersurface in the corresponding ambient space and so, it can be 
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defined by a single polynomial lying in the ideal (F^ a \ G^) (resp. (F^^ 1 ^, G^ " -1 !)) whose 
total degree (see [I0l Lemma 2]) is bounded by degVi (resp. degVo). 

Applying the Bezout Inequality (see \10\ Theorem 1]) to obtain an upper bound for 
deg Vi (resp. degVo), we conclude that there exist polynomials Pi, for i = 1, . . . , s (resp 
i = s + 1, . . . , n + r) meeting the conditions of Corollary [3D whose total degrees can be 
bounded by degPj < <i( CT + 1 )( n + r ) (resp. deg-Pj < <i°"( n+T ')), where d denotes an upper bound 
for the total degree of the polynomials in the input system. 

A similar result can be obtained in the case of Theorem [25] by considering the linear 
projections 

whose domains are the irreducible varieties Vi and Vo defined in (I16p . respectively. 

Analogously, the projections which allow us to estimate the degrees of the polynomials 
in Theorem 1281 are 

for i = 1, . . . , r, all of them defined over the irreducible variety 

Vi := = 0}. (17) 

5.3.2 Algorithmic Issues 

This section presents algorithmic procedures for the computation of the following objects: 

(1) The differentiation index of the system ([1]). 

(2) A differential transcendence basis of a first-order system of type (pQ) preserving the 
order. 

(3) The implicit relations in separated variables given by Theorems [25] & [28] and Corol- 
lary EH 

For algorithmic reasons, we will assume that our base field k is either Q or Q(i) with 
the usual derivation (in the case when k = Q(t), the input polynomials will be assumed 
to have coefficients in Q[i]). 

We start with a brief description of the computational model. 

Basic algorithmic notions 

The objects our algorithms deal with are multivariate polynomials which will be encoded 
by means of straight-line programs (i.e., arithmetic circuits which enable us to evaluate 
them at any given point). The number of instructions in the program is called the length 
of the straight-line program. For a brief description of the algorithmic model and the data 
structure we will use, we refer the reader to Section 2.2] and the references therein. 

The basic subroutine we use is a polynomial-time probabilistic procedure for the com- 
putation of the rank of a matrix with polynomial entries. Roughly speaking, the problem 
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is reduced to the computation of the rank of a matrix with entries in Q by randomly 
choosing integer values for the variables and evaluating all matrix entries at them (see 
Lemma 24] for the error probability analysis). 

Our algorithms take as input a straight-line program of length L encoding the poly- 
nomials /x,...,/ n , g±,...,g r appearing in the system (pQ). However, the intermediate 
computations involve not only these polynomials but also their successive derivatives and 
so, straight-line programs for these derivatives are needed as well. The existence of "short" 
straight-line programs encoding them is ensured by [7J Lemma 21] (see also |21( Section 
5.2]). 

We recall that all our computations should be performed over the differential field K 
with the derivation induced by Xi = fi {i = 1, . . . , r) but, since the number of derivatives 
involved in each computation is controlled, they can be achieved over the polynomial rings 
k[X, U™] for adequate choices of I. We assume that an upper bound d € N for the degrees 
of the polynomials f±, . . . , f n , gi, . . . ,g r is known. In order to estimate complexities and 
error probabilities, we also need upper bounds for the degrees of the polynomials obtained 
by successive differentiation of the input polynomials (i.e., those giving the isomorphism 
between Frac(/c{Y, X, U}/A) and K), which can be found in Notation 6 and Remark 
25]. 

Computation of the differentiation index 

According to Definition [TS] and Notation [T31 we have a := min{&; 6 No / ft = Mfc+i}- 
This minimum is obtained by computing and comparing the ranks of the matrices Zk,e-i 
(which are computed over polynomial rings) for successive values of k 6 N. The algorithm 
finishes, since we have an a priori upper bound for a (see ©). 

The previous computation of a can be achieved with error probability bounded by e 
within a complexity which is polynomial in n,m,r, and linear in logd, loge and L (recall 
that d and L are upper bounds for the degrees and the size of the straight-line program 
representation of the input polynomials). Note that this algorithm can also be applied 
for the computation of the modified index a introduced in Subsection 15.11 with the same 
complexity bounds. 

Computation of a differential transcendence basis for first-order systems 

The algorithmic computation of a differential transcendence basis of the differential field 
extension L <^-> Frac(L{X, U}/A) preserving the order after localization follows the pro- 
cedure underlying the proof of Proposition [30t 

• compute a transcendence basis Bo C {X±, . . . , X n , U±, . . . , U m } of L ■— > !Fq\ 

• choose a subset Bo C Bo with m — r elements such that B\ := Bo U Bo is a transcen- 
dence basis of L ■— > T\ . 

Then, the variables W\, . . . , W m - r whose derivatives lie in Bo form a differential transcen- 
dence basis of L > Frac(L{X, U}/A) with the required property. 

The set Bo is constructed recursively by adding one variable at a time. In order to 
determine whether a subset of variables in To (resp. T\) is transcendental over the field L, 
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we use the fact that A a j A a (resp. T\ A a+ \ / A. a+ \) . Thus, the problem amounts 

to determine whether a subset of variables in a quotient of a polynomial ring by a prime 
ideal is transcendental over the base field, which is done by applying the Jacobian criterion 
from commutative algebra (see Lemma 19]). 

Computation of the implicit equations 

As we have shown in Subsection 15.3.11 the polynomials Pi's of Theorems [25] & [28] and 
Corollary [31] can be interpreted as eliminating polynomials of the image of the algebraic 
varieties defined in (|16[) and (|17|) under suitable linear projections. Therefore, they can be 
computed by means of an algorithm based on standard algebraic elimination procedures 
(see [11] and [32]). For simplicity, we assume k := Q. The following complexity result can 
be obtained: 

Proposition 33 There is a probabilistic algorithm which computes the polynomials Pi of 
Theorems{25\ & \2^ and Corollary \31\ with error probability bounded by e, with < e < 1, 
and within complexity 0(log(l/e)d 2 L) II(n + m, maxj deg Vj)), where U is a suitable two- 
variate universal polynomial. ■ 

We omit the proof of this result in the present article, since it is rather long and 
technical, and follows closely the proof of [7] Proposition 46]. 
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